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Abstract 

We investigate the ground-state energy and spin of disordered quantum 
dots using spin-density-functional theory. Fluctuations of addition energies 
(Coulomb-blockade peak spacings) do not scale with average addition en- 
ergy but remain proportional to level spacing. With increasing interaction 
strength, the even-odd alternation of addition energies disappears, and the 
probability of non-minimal spin increases, but never exceeds 50%. Within a 
two-orbital model, we show that the off-diagonal Coulomb matrix elements 

help stabilize a ground state of minimal spin. 
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The control of spin in semiconductor nanostructures is essential for a number of 
applications such as spintronics |2| and quantum bits and for fundamental studies of the 
Kondo effect in quantum dots [|]]. In clean quantum dots with circular symmetry, Hund's 
rule is observed - high-spin states (i.e. states of non-minimal spin) appear for partly filled 
shells of degenerate single-particle levels [[J . In disordered or chaotic quantum dots , high- 
spin states are suppressed by the rarity of degenerate or nearly degenerate levels. However, 
the observed absence of even-odd alternation of addition energies in quantum dots [f^-|IOH 
is consistent with ground states of non-minimal spin. Recently, evidence for such high- 
spin states has been obtained from studies of the magnetic dispersion of Coulomb-blockade 



peaks |TT,12|. A number of theoretical treatments have considered the conditions under 



which high-spin ground states might occur [|13| -|2"Ci|j . Many of these works [|16|.|T^.[nj address 



the regime of large dimensionless conductance g » 1 ||21||: other works rely on the validity 



of the Hartree-Fock approximation p^ , [l7| . Simulations on small lattice models indicate a 
significant fraction of S = 1 ground states flT5| , but we infer g ~ 0.1 — 0.3 in the simulations, 
compared to the experimental range g > 1 pHID . 



Density-functional theory enables in silico experiments on the ground-state spin of small, 
disordered quantum dots for which g > 1, over a wide range of interaction strengths. The 
results presented here challenge the view that high-spin states dominate in disordered dots 



for r s > 1 [IjJ: even-odd alternation of addition energies disappears by r s = 1.25, but the 
probability of a high-spin ground state never exceeds 50%. Using a two-orbital model, we 
show that off-diagonal Coulomb matrix elements help stabilize a ground states of minimum 



spin, as argued by Jacquod and Stone p0| . We also find that addition-energy fluctuations 



do not scale with the average addition energy, but remain proportional to the single-particle 



level spacing, consistent with our previous results for spin-polarized dots |25 |. 

The ground-state energy and spin of disordered, two-dimensional quantum dots are ob- 
tained within spin-density-functional theory (SDFT). Specifically, we solve the following 



Kohn-Sham equations numerically, and iterate until self-consistent solutions are ob- 
tained; 



2m* k 7 |r — r'| 8pP{r) 2 

where the density is p(r) = So-P CT ( r ) — So- Si l^( r )| 2 > an d the sum is taken over the N 
lowest-energy orbitals. Here a denotes the spin index, £(r) is the local spin polarization, 
and E xc [p,(] is the exchange-correlation energy functional p3| . We use the effective mass 
for GaAs, m* = 0.067m, and a two-dimensional harmonic-oscillator confining potential with 
huo = 3.0meV. The dimensionless interaction strength is measured by (e 2 / K£ )/hu: Q or r s 



2% and is controlled by changing the dielectric k, where k = 12.9 for GaAs. 

We use the same parameters for the total impurity potential VJ mp (r) = 
X^j(7i/27rA 2 )exp (— |r — rj| 2 /2A 2 ) as for a previous study of spin-polarized dots [P5| . The 
density of individual impurities is nj mp = 1.03 x 10 _3 nm -2 ; each has a Gaussian poten- 
tial profile, with strength j t distributed on [— W/2,W/2] with W = 10h 2 /m*, and width 



A = 4/(2^2), where £ = \Jh/w*uj Q ~ 19.5 nm is the quantum-confinement width of the 
parabolic potential. The resulting mean free path, I = vpT ~ 120nm, is comparable to the 
dot diameter L = 120 — 160nm, where L increases with interaction strength. Thus the dots 



are marginally in the ballistic regime and have a dimensionless conductance g ~ 2 [21 



The density functional method is known to give accurate ground-state energies and 



spins for a clean dot [^(J . Comparison with quantum- Monte-Carlo calculations [ 27 1 confirms 
that SDFT is valid for interaction strengths up to (e 2 / kIq) / hujQ = 6 (r s ~ 8) and up to 8 
electrons. For disordered dots, we have confirmed the accuracy of the SDFT results by exact 
diagonalization studies of three electrons, over the range < (e 2 / k£ ) / tku = 2.86 (0 < r s < 
2.92). 

At low temperatures, electron hopping into a dot is suppressed except when the ground- 
state free energies for iV — 1 and N electrons are degenerate. This condition determines the 
position of the iVth conductance peak. The iVth peak spacing or addition energy is given 
by A(JV) = E(N + 1) -2E(N) + E(N -1), where the ground-state energy E(N) is obtained 
from 
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W = E<"£/f^-E/™^* + S«. (2) 

Figure 1 shows (a) the average addition energy (A) and (b) its fluctuation 5 A = ((A — 
(A)) 2 ) 1 / 2 , as functions of electron-electron interaction strength (e 2 / k£q) / Tlujq for AT = 10 and 
N = 11. Dashed curves are for total spin fixed at its minimum, S = or 1/2. Solid curves 
are for unrestricted spin, i.e. the lowest energy E(N) is found among all possible spin values. 
Also shown in (a) is the average noninteracting level spacing ( A ) and in (b) its fluctuations 
(SAq), for dots of the same size as in the interacting case. In Fig. 1(a), there is an even- 
odd alternation in addition energies in the weak-interaction regime, (e 2 / ' kIq) /Tiujq < 0.95 
(r s < 0.76), which reflects the energy cost of adding every odd electron to a new orbital 
level. No alternation is observed in the regime of strong interactions (r s > 1) where (A) 
increases nearly linearly function of interaction. 

In Fig. 1(b), the magnitudes of the addition-energy fluctuations for N = 10 and AT = 11 
merge around r s ~ 1.25. The interacting fluctuations 5A are always smaller than the 
fluctuations 5Aq of the noninteracting level spacing, where the latter satisfy the random- 
matrix-theory relation for the Gaussian orthogonal ensemble |28f 5A ~ — 1(A ) — 
0.52(A ). The interacting fluctuations are roughly 5 A m 0.28(A ), and clearly do not 
scale with the average addition energy (A). In fact, the interacting fluctuations are slightly 
smaller than those for the same number of interacting spin- polarized electrons in the same 
potential |2~5| . (In the spin-polarized case, the addition-energy fluctuations were shown to 



be dominated by the single-particle level-spacing fluctuations because of the effectiveness of 
screening in reducing the Coulomb fluctuations.) 

Comparing Figs. 1 (a) and (b), one sees that in the range r s = 0.76 — 1.25, the average 
addition energy is the same for AT = 10 and A" = 11, but the A^ = 10 distribution has a 
significantly larger width. For r s = 1.25 and below, the shapes of the A^ = 10 and AT = 11 



distributions are qualitatively different |29] consistent with recent experiment observations 



for a high-electron-density GaAs quantum dot with r s = 0.72 UTT] . 



For the strong- interaction regime r s > 1, fluctuations in addition energy are reduced 



when the spin states are unrestricted. In the lower right inset, we show the distributions of 
A(iV) for N = 10 (upper) and for N — 11 (lower) for interaction strength (e 2 / k£q) /Tiujq = 
1.91(r s = 1.77), both for restricted and unrestricted spin. The distributions of A are roughly 
Gaussian, which agrees with experiments at r s > 1 [pHIO]- (For comparison, the Wigner- 



Dyson distribution [^] for the noninteracting case is shown in the upper left inset.) The 
distributions of A become narrower in the spin unrestricted case: the low-energy tail of 
the N = 10 distribution is reduced, while the high-energy tail is reduced for N = 11. 
We trace this behavior to the appearance of S = 1 ground states for 10 electrons (see 
below). Since A(N = 10) = E(ll) — 2E(10) + E(9), the appearance of new, low values 
for E(10) removes some of the lowest values of A(N = 10). Similarly, since A(iV = 11) = 
£"(12) — 2£'(11) + -E'(IO), new, low values of E(10) remove some of the highest values of 
A(N = 11). 

Figure 2(a) shows the probabilities of the different ground-state spins S versus electron- 
electron interaction strength. Solid curves are for N = 10 electrons (integer spin) and dashed 
curves are for N = 11 electrons (half-integer spin). With increasing interaction strength, the 
probabilities of ground states of S = 1, 3/2, and 2 increase, while the probabilities of even 
higher spin states (S = 5/2, etc.) remain negligible. The probability of S = 1 for iV = 10 is 
always higher than that of S = 3/2 for N = 11. The probability of an S — 1 ground state 
reaches a maximum around (e 2 / k£q,) /Uujq ~ 1.91(r s ~ 1.77), and never exceeds 50%. In the 
inset, we show the increase with interaction strength of the probability of spin blockade, i.e. 
subsequent ground states with \ AS\ > 1/2, which leads to a suppressed Coulomb-blockade 
peak |0|. 

Diagonalization studies on small lattices also find a significant fraction of S = 1 ground 
states, with a smaller likelihood of S = 3/2 [[0J. The onset of high-spin ground states 
occurs at smaller r s ~ 0.2, and addition-energy fluctuations are larger, scaling as 10-20% of 
the average addition energy. We attribute these differences to weak screening in the lattic 
models due to the low dimensionless conductance of g = 0.1 — 0.3, compared to g ~ 2 in our 
dots. 



Berkovits (T^j has argued that high-spin ground states are favored not only by exchange 
energy (which favors spin alignment) but also by the enhanced Coulomb repulsion between 
two electrons in the same spatial orbital. Opposing these effects is the single-particle energy 
cost of promoting an electron to a new orbital. He observes that it is much more likely 
to find two orbitals close in energy, producing an S = 1 ground state, than to find three 
orbitals close in energy, as required for an S = 3/2 ground state. 

These arguments are consistent with our SDFT results up to r s ~ 1, but do not account 
for the observed saturation of the probability of high-spin ground states at larger r s . To 
understand this saturation, we consider the possible states of two electrons occupying two 
spin-degenerate orbitals near the Fermi energy. (The other electrons in the dot are assumed 
to pair up in lower-energy orbitals.) For this two-orbital model, there are three degenerate 
S = 1 states consisting of one electron in each of the two orbitals, and three nondegenerate 
S = states. 

To evaluate the energies of the competing S = and S — 1 two-electron states, we use 
the basis of single-particle eigenstates 0°(r) and 0° +1 (r) with energies e n and e n+ i of non- 
interacting electrons. The various diagonal and off-diagonal Coulomb matrix elements are 
evaluated in the random-phase approximation (RPA) which accounts for the screening effect 
of the other electrons in the dot. The energy of the three degenerate S = 1 states is E(S = 
1) = e n + e n+1 + U nin+1 -X n>n+1 , where U n>n > = e/^ n)ri (r)p°, n ,(r)dr is the screened Coulomb 
interaction between two electrons in orbitals n and n', and X n ^ n > = e J (p n y{r)p^, n {v)dr is 
the screened exchange interaction. Here, <f n ,n'( r ) is the screened potential due to an electron, 
which is evaluated in Fourier representation as <p n ,n'(<l) = (2vre//t|q|)(p° n /(q)/e(q)), where 
Pnn'( r ) = 0n( r )^n'( r )- The dielectric function e(q) is approximated as e RPA (q) = l — v q x{q) 
where v q = 2n/q and the susceptibility x(g) = — (e 2 / K>)(dn/dfi)F(q/2kF). The Lindhard 
polarizability for 2D is F(x) = 1 for x < 1 and F(x) = 1 — \Jl — (1/x) 2 for x > 1. 

The energy E(S = 0) of the lowest 5 = state is obtained by diagonalizing the following 
3x3 matrix; 
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Hs=o 



2Cn ~\~ U ni n v2^n,n,n,n+l -^n,n+l 

,n+l + X ntn+ i \^2U n+ i 

,n+l,n+l,n 

l,n+l,n+l,ra 2e n+ i + Z7 n +l,n+l 



(3) 



where the off-diagonal Coulomb matrix elements are U n ,n,n,n' — e/ ^n,n(r)p° n'( r )^ r - We 
find that the magnitudes of the U n ,n,n,n' are comparable to the exchange energy X n n+1 
|20| . It is seen in Fig. 2(b) that for the two-orbital model the average of AE = E(S = 
0) — E(S = 1) agrees reasonably well with our SDFT results for all strengths of interaction. 
In contrast, placing the two electrons in the lowest single-particle orbital 0°(r), gives an 
energy 2e n + U nn which is significantly larger than E(S = 0) at large r s . Even the Hartree- 
Fock (HF) approximation for i.e. the best doubly-occupied single orbital, overestimates 
the energy of the lowest 5 = state and fails to account for the saturation of the S = 1 
probability at large r s . It is evident that for the two-orbital model the off-diagonal Coulomb 
matrix elements help stabilize the S = ground state by creating a low-energy hybridization 



of the three 5 = basis states, as proposed by Jacquod and Stone |20| . 

We acknowledge fruitful discussions with B. Altshuler, R. Berkovits, D. Goldhaber- 
Gordon, B. Halperin, C. Marcus, M. Stopa, and F. Zhou. 
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Figure 1: (a) Average addition energy (A) and (b) fluctuation 5 A = ((A — (A)) 2 ) 1 / 2 , as functions 
of electron-electron interaction strength (e 2 / k£ ) /hcu for N = 10 and N = 11. Solid 
curves are for unrestricted spin; dashed curves are for total spin fixed at S — or 1/2. 
For each data point, the disorder average is taken over more than 1000 different impu- 
rity configurations. The average noninteracting level spacing and its fluctuations are 
shown for the same size dots as in the interacting case. Lower right inset - Distribu- 
tion of addition energies A for N = 10 (upper) and for N = 11 (lower) for interaction 
strength (e 2 / k£ ) / 'Tuvq = 1.94. Upper left inset - Distribution of noninteracting addi- 
tion energies (level spacings) for dots of the same size as in the lower inset. 

Figure 2: (a) Probability of a spin S ground state as a function of electron-electron interaction 
strength (e 2 / K,£ )/huj . Solid curves are for iV = 10 electrons (integer spin) and dashed 
curves are for N = 11 electrons (half-integer spin). The circles give the probability 
of a spin 5 = 1 ground state for the full two-orbital model. Inset - Probability of 
spin blockade, i.e. subsequent ground states with \AS\ > 1/2. (b) Average total 
energy difference AE between states with S = and S = 1 as a function of electron- 
electron interaction strength (e 2 / ' k£q) /Huq. The solid curve shows the SDFT results. 
Also shown are results of the two-orbital model: exact (lowest dashed curve), Hartree- 
Fock (central dashed curve), doubly-occupied lowest orbital (highest dashed curve), 
with parameters evaluated for the 5th and 6th non-interacting orbitals. The apparent 
agreement between the HF result and SDFT at r s ~ 1 is not meaningful - the HF 
result is only an approximation to the exact result (lowest dashed curve). Inset - 
Average level spacing and screened Coulomb interactions for two-orbital model. 
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